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Abstract 



In a previous paper (Corrigan-Sasaki), many remarkable properties of classical 
Calogero and Sutherland systems at equilibrium are reported. For example, the mini- 
mum energies, frequencies of small oscillations and the eigenvalues of Lax pair matri- 
ces at equilibrium are all "integer valued" . The equilibrium positions of Calogero and 
Sutherland systems for the classical root systems (A r , B r , C r and D r ) correspond to 
the zeros of Hermite, Laguerre, Jacobi and Chebyshev polynomials. Here we define and 
derive the corresponding polynomials for the exceptional (Eg, E-j, Eg, and G%) and 
non-crystallographic (^(m), #3 and H4) root systems. They do not have orthogonality 
but share many other properties with the above mentioned classical polynomials. 



1 Introduction 



The relationship between classical and quantum integrability has fascinated many physicists 
and mathematician. In a recent paper by Corrigan and Sasaki HJ, this issue has been 
extensively investigated in the framework of Calogero-Moser systems @, 0, §]• One major 
result is that certain "quantised" information seems to be encoded in the classical system. 
For example, the eigenvalues of classical Lax pair matrices at the equilibrium points are 
"integer valued" [Qj • The connection between the zeros of Hermite and Laguerre polynomials 
and the equilibrium points of A r and B r (D r ) Calogero systems has been known for many 
years H H 0- In it is found out that the zeros of Jacobi polynomials are related to the 
equilibrium points of BC r (D r ) Sutherland system. In the present paper we define and derive 
the polynomials associated with the equilibrium points of the other Calogero and Sutherland 
systems. These are associated with Calogero systems based on non-crystallographic root 
systems, Calogero and Sutherland systems based on the exceptional root systems and the A r 
Sutherland systems. The Chebyshev polynomials ( [5 .3D are associated with the A r Sutherland 
systems. 

In general, the polynomials are determined by the potential, q 2 + 1/q 2 (the Calogero 
system 0) and l/sin 2 g (the Sutherland system 0), the root system A and the set of 
weights TZ. For the classical root systems and for the (non-trivial) smallest dimensional TZ, 
that is the set of vector weights V or the set of short roots A$, the polynomials turn out to 
be classical orthogonal polynomials; Hermite, Laguerre, Jacobi and Chebyshev polynomials 
Pj. The orthogonality does not hold for the polynomials for the exceptional root systems 
and for the classical root systems with generic TZ's. Like their classical counterparts, these 
new polynomials have "integer coefficients" only, if multiplied by a certain factor. In most 
cases, it is possible to define the polynomials to be monic (that is, the highest degree term 
has unit coefficient) and integer coefficients only. Some polynomials are too lengthy to be 
displayed in the paper; an E% polynomial has 121 terms and its typical integer coefficient 
has about 150 digits. They are presented in ||. Some root systems are related by Dynkin 
diagram foldings; A 2r -i — > C r , D r+ i — > B r , E e — > F 4 and D4 — ► G 2 - These imply relations 
among the corresponding Calogero-Moser systems at certain ratio of the coupling constants. 
These, in turn, imply relations among the corresponding polynomials, which are determined 
independently. These relations are either identities among classical polynomials, many of 
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which are "new" in the sense they are not listed in standard mathematical textbooks ||, 
or they provide non-trivial checks for the newly derived polynomials. The significance and 
other detailed properties of these new polynomials deserve further study. 

This paper is organised as follows. In section two a brief introduction of Calogero- 
Moser systems is given to set the stage and notation. Equations for determining equilibrium 
points are discussed in some detail. In section three Coxeter (Weyl) invariant polynomials 
associated with equilibrium positions are introduced for a set of weights 1Z for Calogero 
and Sutherland systems. For the rational potential (Calogero systems) the definition is 
almost unique, whereas we have several choices of the definitions of the polynomials for the 
trigonometric potential (Sutherland system). Section four and five are the main body of 
the paper. The Coxeter (Weyl) invariant polynomials are determined and presented for all 
root systems A and for major choices of 7£'s for Calogero (section four) and Sutherland 
systems (section five). Section six is for summary and comments. We will present a heuristic 
argument for deriving the classical orthogonal polynomials starting from the pre-potentials 
(|2.4j) of Calogero and Sutherland systems. 



2 Equilibrium in Calogero- Moser System 

Let us start with a brief introduction of Calogero-Moser systems 0, |3|, Hj. We stick to the 
notation of a recent paper |]], unless otherwise mentioned. Calogero-Moser systems are in- 
tegrable multiparticle dynamical systems at the classical as well as quantum levels. They 
have a long range potential (rational, trigonometric, hyperbolic and elliptic) and the inte- 
grable multiparticle interactions are governed by the root systems ||10|| . Classical integrability 
through Lax formalism is known for all potentials for the classical root systems [|1(J as well 



as for the exceptional [II], |T^[ and non-crystallographic [TJ|] root systems. Quantum integra- 
bility of the systems having degenerate potentials (rational, trigonometric and hyperbolic) 
is now systematically understood for all root systems in terms of Dunkl operator formalism 



13 , 14 1 and the quantum Lax pair formalism [15, [If]]. To a system of r particles in one 
dimension, we associate a root system A of rank r. This is a set of vectors in W invariant 
under reflections in the hyperplane perpendicular to each vector in A: 



2a 

A 3 s a {(3) = (3 - (a v ■ /3)a, a v = — , a,^A. (2.1) 

ar 
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The set of reflections {s a \ a G A} generates a finite reflection group Ga, known as a Cox- 
eter (or Weyl) group. Among Calogero-Moser systems the Calogero systems (with q 2 + 1/q 2 
potential) and the Sutherland systems (with 1 / sin 2 g potential) have discrete energy eigenval- 
ues only when quantised. The Calogero and Sutherland systems have equilibrium positions, 
which are characterised by two equivalent ways [l] . That is where the classical potential takes 
the absolute minimum and simultaneously the groundstate wavefunction takes the absolute 
maximum. At the equilibrium positions of the Calogero and Sutherland systems, associated 



spin exchange models are defined for each root system |fL7|j . The best known example is the 
Haldane-Shastry model which is based on A r Sutherland systems [JT3]. The integrability and 



the well ordered spectrum of the spin exchange models are closely related with the special 
properties of systems at equilibrium |1J. 

The classical Hamiltonians of the Calogero and Sutherland systems readQ 



n 



c 



:P 2 + V c , V t 



u 2 „ 1 



9 y 



2 ^ (p-q) 2 



c 



9 < d 



9 y 



(2.2) 



2 ^ + sin (p-q) 



In these formulae, A + is the set of positive roots and uo > is the angular frequency of the 
confining harmonic potential, g p > are real coupling constants which are defined on orbits 
of the corresponding Coxeter group, i.e., they are identical for roots in the same orbit. The 



classical potential Vc can be written succinctly in terms of a pre-potential W [15 



V 



c 



1 y> fdWV 



(2.3) 



in which 



W 



UJ 



P eA + 
g P ^g\sin(p- q)\, 

P eA+ 



(2.4) 



and So is the minimum energy. Let us recall that the pre-potential W is related to the ground 



e w (eq.(2.6) of 



state wavefunction of the quantum theory (fio by (fio = 

1 For A = BC r the trigonometric potential should read gf d Y^ p eA M ^ I s ^ (p'l)~' r 2 9l S p g a 
?) + 9s(gs + 2g L )/2 E P eA s+ V sin2 (-° ' ?)« with At = 2 > A = 4 and Ps = L 



and that W, V c 

2 ^ l/sin 2 (p- 
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and Tic are Coxeter (Weyl) invariant: 

H c (p,q)=Hc(s a (p),s a (q)), W(q)=W(s a (q)), V c (q) = V c (s a (q)) (Va G A). (2.5) 
The classical equilibrium point 

p = 0, g = g (2.6) 



is determined by the equations JT 



or equivalently 



= (j = l,...,r). (2.7) 
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(j,Z = l,...,r), (2i 



In other words, it is a minimal point of the classical potential Va, and simultaneously it is a 

maximal point of the pre-potential and of the ground state wavefunction <f) = e w , since 

the matrix determining the frequencies of small oscillations around the equilibrium 

d 2 W 
dqjdqi 

is negative definite [lj. The equilibrium points are not unique. There is one equilibrium 
point in each Weyl chamber (alcove) |1], that is if q is an equilibrium point, so is s p (q), 
Vp G A, due to the Coxeter (Weyl) invariance of W (|2.5|) . It is also easy to see that if q is 
an equilibrium point, so is —q. 

The equilibrium equation for the pre-potential W, for Calogero systems based on simply 
laced root systems, that is A r , D r , E r , I 2 (odd) and H r , reads: 

P v - 

— = -i- 

P eA+ 

If we define a rescaled equilibrium point by 



v P ' q 9 

PGA+ y 



q=J-q, (2.9) 



it satisfies a simple equation independent of the coupling constant: 

P 



P6A+ * 

For Calogero systems based on non-simply laced root systems, that is B r , C r , G2 and 
/2(even)0, the equation reads: 



2 For ^(even) we have k = g e /go- 



Again a rescaled equilibrium point 

q= x [^q, (2.11) 
V 9l 

satisfies a simple equation depending only on the ratio of the two coupling constants gs and 
9l- 

£ - £ -- + k £ A: = 9, iE* (2.12) 



As is clear from ( 2.1U| ) and ( [2.12 ), the equilibrium point q (q) is independent of the normal- 



isation of roots in A. 

The situation is simpler in the Sutherland systems which do not have an extra parameter 
uj. For crystallographic simply laced root systems, that is A r , D r and E r , the equation for q 
is independent of the coupling constant: 

^pcot(p-g) =0. (2.13) 

P6A+ 

For crystallographic non-simply laced root systems, that is B r , C r , F4 and G 2 , the equation 
for q depends only on the ratio of the two coupling constants gs and g^- 

pcot(p- q) + k P cot (p- q) = ° ; k=— . (2.14) 

For the BC r system, which has three coupling constants gs, gu and gi for the short, middle 
and long roots, the equation depends on two coupling ratios: 



Y. pcot(p • q) + ki V" pcot(p • q) + k 2 V" p cot(p • q) = 0, k x = — , k 2 = — . 

P6A M+ P eA 3+ P eA L+ 9M gu 

(2.15) 



3 Polynomials 

Here we give the general definitions of the Coxeter (Weyl) invariant polynomials associated 
with equilibrium positions in Calogero and Sutherland systems. Naturally, the definitions 
for the Calogero systems are different from those for the Sutherland systems except for the 
common features that the polynomials are Coxeter (Weyl) invariant and are specified by the 
root system A and a set of D vectors 71 

71 = {p (1) ,...,/i (D) |/i (a) GR r }, (3.1) 
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which form a single orbit of the corresponding reflection (Weyl) group The set of values 
at the equilibrium, {// • q\ fi 6 7Z}, is Coxeter (Weyl) invariant. In this paper we consider 
only such TZ's that are customarily used for Lax pairs. They are the set of roots A itself 
for simply laced root systems, the set of long (short, middle) roots A^ (As, Am) for non- 
simp ly laced root systems and the so-called sets of minimal weights. The latter is better 
specified by the corresponding fundamental representations, which are all the fundamental 
representations of A r , the vector (V), spinor (S) and conjugate spinor (S) representations 
of D r and 27 (27) of E 6 and 56 of E 7 . 

For Calogero systems the definition is rather unique and it is given by 

P%(k\x) = Y[(x-(i-q), (3.2) 

in which k denotes the possible dependence on the ratio of the coupling constants, for 
the systems based non-simply laced root systems Q2.12 ). It should be noted that the above 
polynomial depends on the normalisation of the vectors fx G 7Z implicitly. Changing 7Z — > cTZ 
(fx — > cfx) can be absorbed by rescaling of x: 

Pf(k\x) = Y[ (x - cfi ■ q) = c D P%(k\x/c). (3.3) 
veil 

For Sutherland systems we have several candidates for polynomials: 

PA,.&\x) = ]l(x-8m(ji-q) > ), Pl s2 (k\x)= JJ (x - sin(2/z • q)) , (3.4) 
fien lien 

Plc(k\x)= Y[(x-cos(fi-q)^, Pl c2 (k\x) = H(x-co S (2fx-q)}, (3.5) 

in which k denotes possible dependence on the ratio(s) of coupling constants, as before. Not 
all of them give interesting objects, as we will see presently. In all cases the polynomials are 
monic and of degree D. 

In case TZ is even, that is, 

fx e JZ -fx e 71, (3.6) 

then sometimes it is advantageous to consider P]±(k\x), P™ s (k\x) and Pa, S 2(^\ x ) as polyno- 
mials in y = x 2 of degree D/2: 

II (y - (a* • <?) 2 ) » U(y- sin2 ^ ■ ?)). U{v- sin2 ( 2 ^ ■ ?))> ( 3 - 7 ) 

IJ.&1+ n€H+ AiG7?.+ 



in which 1Z + is the positive part of 7Z. In this case the "cosine" polynomials P& c r 2 )(k\ 



x 



should better be redefined as 

Pl!Mx) = J] (x - cosOi • qj), PlUk\x) = J] (x - cos(2/i ■ g)), (3.8) 



since the original polynomials ( |3.5| ) are the squares of the new ones. It is easy to see that 
P]± a (k\y) and P^ + c2 (k\x) are equivalent: 

F%.(k\x) = (-2)-^ 2 Pj: 2 (fc|l - 2a; 2 ). (3.9) 

Likewise, for even TZ, P£ s2 {k\x) is a "square" of P^ + c2 (k\x): 

P*Mx) = J] ( x - sin ( 2 / x • <?)) = II (^ 2 -sin 2 (2/i-g)^ 

v /v — cos(2/i • g) j , »" = !—./• 



] J ru — cos(2/i • g)J (^— u — cos(2/i • q) 

p2; C 2(k\u)p2; c2 (k\-u). (3.10) 



The right hand side is an even polynomial in u, thus it is a polynomial in u 2 and in x 2 . 
The change of variables u <-> x corresponds to the change in the character of the variables, 
cos <-» sin. This imposes a quite non-trivial check for the s2 and c2 polynomials which are 
determined separately. 

As shown in the following sections, the polynomials associated with the classical root 
systems (A r , B r , C r and D r ) and ^(tu) are either classical polynomials for the smallest 
dimensional 71 or those closely related to them, see for example, ( f4.32|) , (|4.33|) , ( |5.41|) , (|5.42|) . 



For the exceptional and non-crystallographic root systems, the equilibrium positions are 
evaluated numerically and the polynomials are obtained by rationalisation of the coefficients 
in terms of Mathematica. At each step, the result is verified by many consistency checks; the 
"integer eigenvalues" of the matrix ( |2.8j ) for the values of q, the identities implied by Dynkin 



diagram foldings and identities ( 3.10 ) for the polynomials. Let us conclude this section with 
an important remark that these polynomials are independent of the specific representation of 
the root and weight vectors. In other words, the polynomials are Coxeter (Weyl) invariant. 

4 Calogero Systems 

Let us first discuss the systems based on the classical root systems. 
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4.1 A r 

The equations ( |2.10| ) for A = A r read 



r+l 



Err = ^ (j = i,...,r + i). (4.i; 

^7 9j ~ * 



Mi 

These determine {g-,- = ^J^Qj \j = 1, • • • , t + 1} to be the zeros of the Hermite polynomial 
H r +i(x) [U, with the Rodrigues' formula 

H n (x) = (-l) n e* 2 e~ x2 = 2 n x n + ■■■. (4.2) 

If ordered by the value, qi > q 2 > ■ ■ ■ > q r +i or reverse, they possess the symmetry 

Qj = -Qr+2-j, (4-3) 
and especially <?( r +2)/2 = for r even. Thus we have 

qi + g 2 + • • • + Qr+i = 0. (4.4) 

4.1.1 K = V for A r 

This case was reported by Calogero a quarter century ago ||. The set of weights of the 
vector representation is 

V = fa = e 3 - -— e« | 3 = 1, • • • , r + l}. (4.5) 

2=1 

Throughout this paper we denote an orthonormal basis of M r (R r+1 for A r case) by {e^}. 

In this case, we have ■ q — q~j due to ( |4.4|) and ji 2 — r/(r + 1). The polynomial Q3.2|) is 
given by the Hermite polynomial 

r+l 

P r v (z) = Pj>) = j[ (x - ^ = 2-^H r+1 (x). (4.6) 

They are orthogonal to each other: 

P r v (x)P s v (x) e" x2 da; oc 5 rs . (4.7) 

Needless to say, Hermite polynomials are of integer coefficients. It is interesting to note that 
another definition 

r+l 

Pg(x) = II (* - 2 «0 = ^+i(*/2) = 2 r+1 PXW2). (4.8) 

i=i 

gives a monic polynomial with all integer coefficients. 



4.1.2 TZ = V; for A r 



The set of weights of the i-th fundamental representation ( i-th rank anti-symmetric tensor 
representation, 1 < i < r) is 



+ • • • + ^ I 1 < h < ■ ■ ■ < ji < r + l}, D = A = r + . M . 



(4.9) 



The above V ( |4.5|) is V = Vi. In this case we have // 2 = i(r + 1 — i)/(r + 1). We can show 
that the polynomial (|3.2| ) 

PYHx)= n (*-(& + ■•■&)) =p£ w - 



(4.10) 



1<j'i<— <ii<r+l 

can be expressed in terms of the coefficients of H r+ i(x) by the same method as given in 
section [4.2.5| , and P 2 A i {x) gives a monic polynomial with integer coefficients. 

Here we report only on V2 because it seems that the other representations (3 < i < r — 2) 
do not provide any interesting results. (For lower rank r, the explicit forms of the polynomials 
Pj i (x) can be found in ||.) Due to ( |4.3j ), eq. ([4.10|) becomes 

= n (*-&+«)) ( 4 - n ) 



l<j<l<r+l 



(r+l)/2 



r/2 



n (* 2 - & - (* 2 - & + 

l<7<i<(r+l)/2 

]> 2 - $ n 2 - & - *) 2 ) (* 2 - & + « 



r : odd 



r : even. 



i=i i<i< / <t/2 

Based on the fact that the zeros of Hermite and Laguerre polynomials are related as seen 
from the formulae ( [4.23 ), this can be expressed by using the polynomials associated with 
the B r Calogero systems in the following way: 



x r P^{l/2\x), PZ(x) = x r - l PXJx)P£(3/2\x 



r-l D V 



(4.12) 



The explicit forms of the functions P B L {k\x) for lower r are given in section 4.2.5 . 
4.1.3 K = A for A r 

We have A = {±(ej — ej) | 1 < j < I < r + 1}, D = r(r + 1) and fi 2 = 2. The polynomial 
has a factorized form: 

z f x 2 - 2 (r = 1) 



l<j<l<r+l 



X 



fr > 2) 



(4.13) 



Another definition P| (x) gives a monic polynomial with integer coefficients. 
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4.2 B r and D r 

Assuming q~j ^ 0, the equations (|2.12|) for A = B r with k = gs/gL read 

They determine {g| = | j = 1, . . . , r}, as the zeros of the associated Laguerre polynomial 
L ( r a) (x), with a = k -1 = g s /g L - 1 i, 1, 0- The Rodri Sues' formula reads 

^M^(£)V*" + «H^ + ---. (4-15) 
For the subcase with g s = 0, that is A = D r , {g| = -^gj | j = 1, . . . , r}, are the zeros of the 



associated Laguerre polynomial || [10 



rL[7 l \x) = -x4-i(^), (4.16) 

for which one of the <]j is zero. This also means that the {gj} of B r for gs/gL = 2 or a = 1 
are the same as the non-vanishing {gj} of D r+ i. This can be understood easily from the 
Dynkin diagram folding D r+ i — > B r . We omit C r case, because C r is obtained from B r by 
interchanging the short (S) and long (L) roots. 

4.2.1 TZ = A s for B r 

Since A s = {±e.,- | j = 1, . . . , r} is even, it is advantageous to consider the polynomials in 
y = x 2 , 0, 

r 

i?*(y) = P£ s (k\x) = l[{x 2 - q]) = (-l) r r\Li a \y), a — k — 1 — g s /g L - 1. (4.17) 

i=i 

They are orthogonal to each other: 

POO 

/ P r As (y)P s As (y) y a e~y dy oc 5 rs . (4.18) 
Jo 

It should be stressed that P As (y), a monic polynomial in y, is also a polynomial in the 
parameter a with a// integer coefficients . 
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4.2.2 K = V for D r 

As in the previous example, V = {±ej | j — 1, . . . , r}, we introduce (y = x 2 , 

r 

P?(y) = Pl(x) = l[{x 2 - %) = {-iyr\L^\y). (4.19) 

3=1 

They are orthogonal to each other: 

P7{y)Pj{y)y- l e-vdy oc / L^\(y)L^\(y)ye-y dy oc 5 rs , (4.20) 
o </o 

in which the identity fl4.16|) is used. Corresponding to the above mentioned Dynkin diagram 

folding D r+ i — > B r and ( |4 . 1 6|) , we obtain 

x 2 P^{2\x) = Pl +l (x) = Pi s +1 (0\x). (4.21) 

4.2.3 A 2r -i — > C r and the relationship between Hermite and Laguerre polyno- 
mials 

As is well-known the Dynkin diagram folding A 2r -i — * C r relates the A 2r _i Calogero system 
to the C r (B r ) system with uo — > 2a;, gsioi) — %9 and giids) — 9-, that is a = —1/2. This 
would imply P^^ar) (PD is equal to P^(\j2\x) (gT7|): 

Pl r _ 1 (x) = P^(l/2\x), (4.22) 

which is equivalent to a well-known formula relating Hermite polynomials and Laguerre 
polynomials (eq(5.6.1) of ||): 

H 2r (x) = (-l) r 2 2 V!4- x / 2 )(x 2 ), # 2r+1 (x) = (-l) r 2 2r+1 r\xL^ 2 \x 2 ). (4.23) 



The former corresponds to k = 1/2 and ( [4.221 ). The latter corresponds to k = 3/2 and 
implies 

Pl r (x) = xP^(3/2\x). (4.24) 

Let us recall the corresponding results in the trigonometric case || |TJ. The polynomial 
Pqc + 02(^1 1 kh\ x ) (|5-20|) is proportional to Jacobi polynomial Pl a,t3 \x) with a = k\ + k 2 — 1 
and (3 — k 2 — 1. For ki = 0, k 2 = 1/2 {k x = 0, k 2 = 3/2) it reduces to the Chebyshev 
polynomial of the first (second) kind. As above, k\ = 0, k 2 = 1/2 corresponds to the 
A 2r -\ — > C r folding. 
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4.2.4 H = S and S for D r 

The spinor S and conjugate spinor S representations of D r are minimal representations with 
D = 2 r_1 and the natural normalisation /i 2 = r/4. For odd r, we have the equality — S = S 
which means Pg (x) = (x) for odd r. In fact, the symmetry of the D r Dynkin diagram 
implies that the same formula holds for even r, too. Here we present P^ (a;) for lower 
members of r: 

Pjj; S ' v (x) = x 2 (-24 + 36x 2 - 12x 4 + x 6 ), (4.25) 
P%f(x) = 25 - 3400a; 2 + 13900a; 4 - 20200a; 6 + 12730a; 8 - 3880a; 10 

+580a; 12 - 40a; 14 + x 16 , (4.26) 
P*f(x) = 2" 16 (951356390625 - 24582413628000a; 2 + 229552540380000a; 4 

-1001859665040000a; 6 + 2271780895320000a; 8 - 2992279237056000a; 10 
+2465846485977600x 12 - 1332743493888000a; 14 + 486926396352000a; 16 
-122431951872000a; 18 + 21351239884800a; 20 - 2577889198080a; 22 
+212745830400x 24 - 11668684800a; 26 + 403046400a; 28 - 7864320a; 30 
+65536x 32 ) . (4.27) 



The equality of the three polynomials for V, S and S in D 4l Q4.25| ) reflects the three fold 
symmetry of the -D4 Dynkin diagram. 

4.2.5 K = A L for B r and D r 

The set of long roots of B r is = {±(ej — ej), ±(ej + e;) | 1 < j < I < r}. The polynomial 
P B L (k\x) can be expressed neatly in terms of the coefficients of the polynomial P^{k\x) 
( 4.17|) . Suppose we have two polynomials in y: 



/ = il(v - ^) = £(-!)'«* v* 1- ', ( 4 - 28 ) 

i=l i=0 

9= n (v-( x i- x i) 2 )(v-( x i+ x jf)- ( 4 - 29 ) 



l<i<j<n 

Let us denote b\ = xf, then we obtain g as a symmetric polynomial in bf 



9= II (l/ 2 -2(6i + 6i)2/ + (6i-6i) 2 ), (4-30) 



1< i<j<n 
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and {<2i} are the basis of the symmetric polynomials in bf. 

a i = Yl h n ■■■!>,■ (4-31) 

l<ji<—<ji<n 

Thus g can be expressed in terms of the coefficients {<2j} of / with integer coefficients. For 
example: 

n = 2 : g = y 2 - 2a x y + a\ - 4a 2 , (4.32) 
n = 3 : g = y 6 - 4aij/ 5 + 2(3a 2 - a 2 )y 4 - 2(2a? - aia 2 - 13a 3 )y 3 

+ (a x + 2a 2 a 2 — 7a 2 , — 2Aaia 3 )y 2 — 2(a 2 — 3a 2 )(aia 2 — 9aa)y 

+ a\a\ - Aa\ - Aafa 3 + 18aia 2 a 3 - 2703. (4.33) 

If / is of rational coefficients, so is g. 

We list Pg r L (k\x) for lower members of r. This includes Pp r {x) as a special case of k — 0. 
As remarked before, they are presented as polynomials in y = x 2 : 

P^(fc|a;) = 4(1 + fc) - 4(1 + fc)y + y 2 , (4.34) 
P^(A;|x) = 108(l + A;)(2 + A;) 2 - 324(l + A;)(2 + A;) 2 y + 9(2 + A;) 2 (41 + 32A;)?/ 2 
-4 (2 + fc) (99 + 88 fc + 16 fc 2 ) y 3 + 6 (2 + fc) (17 + 8 fc) 

-12(2 + fc) ? / 5 + ? / 6 , (4.35) 
P^{k\x) = 27648 (1 + fc) (2 + fc) 2 (3 + fc) 3 - 165888 (1 + fc) (2 + fc) 2 (3 + kf y 
+4608(2 + fc) 2 (3 + fc) 3 (91 + 82 fc) y 1 
-512 (2 + fc) (3 + fc) 3 (2282 + 2777 fc + 792 fc 2 ) y 3 
+ 192 (2 + fc) (3 + fc) 2 (15462 + 20235 fc + 8336 fc 2 + 1088 fc 3 ) y A 
-768 (2 + fc) (3 + fc) 2 (2085 + 2167 fc + 688 fc 2 + 64 fc 3 ) y 5 
+64(3 + fc) 2 (17634 + 22113 fc + 9480 fc 2 + 1536 fc 3 + 64 fc 4 ) y 6 
-768(3 + fc) 2 (342 + 327 fc + 96 fc 2 + 8 fc 3 ) y 7 

+48 (3 + fc) (2514 + 2465 fc + 784 fc 2 + 80 fc 3 ) y 8 (4.36) 
-64 (3 + fc) (186 + 123 fc + 20 fc 2 ) y 9 + 240(3 + fc)V° - 24 (3 + fc) y 11 + y 12 . 

As remarked above, P^{k\x) is a polynomial in y and in fc with all integer coefficients 
and is monic in y. The explicit forms of the polynomials P^ r L (k\x) (r = 5,6) and P^ r {x) 
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(r = 4,5,6) can be found in ||. The Dynkin diagram folding D r+ i — > B r relates the 
polynomials 

p£(2\x) {Pi s {2\x)) 2 = Pi +i (x) = Pi L +1 (0\x), (4.37) 
which is the root version of the identity (|4.21 ). 



Next we discuss the systems based on the exceptional root systems. For these we have 
relied on numerical evaluation of the equilibrium points by Mathematica. Large enough 
digits of precision is maintained in internal computations, e.g., we keep 2048 digits for E s 
Sutherland system. We have verified in each case that the fit of the polynomial with rational 
coefficients gives no detectable errors within the working precision. 

4.3 E r 

The E series of the root systems, Eq, Ej and E% are simply laced. The corresponding 
polynomials do not contain any coupling constants. 

4.3.1 K = 27 and A for E 6 

Polynomials for 27 and A, 

pgw = n c* - /* • 9) ^ 2 = 4 a p 2 = 2 )' ( 4 - 3g ) 

peA 

are slightly simplified for a different normalisation of /i e 1Z: 

y/T/327 



e 6 (*) = 3- 27 / 2 P|J(V3x) = J] (x - jk ■ q) (ft = fi/y/3, fi 2 = 4/9) 

= x 3 (200 - 3600 x 2 + 24600 x 4 - 83980 x 6 + 162945 x s - 192840 x 10 

+144876 x 12 - 70416 x u + 22170 x 16 - 4440 x 18 + 540 x 20 - 36 x 22 + x 24 ) , (4.40) 

V3A 



Pg ' {*) = 3- 36 P# 6 (V3x) = H(x-p-q) ( P = P /V3, p 2 = 2/3) 

81920 - 1474560 x 2 + 8970240 x 4 - 22749184 x 6 + 28505088 x 8 - 19829760 x 10 



+8239872 x 12 - 2128896 x 14 + 346944 x 16 - 35328 x 18 + 2160 x 20 - 72 x 22 + x 24 
x (200 - 3600 x 2 + 24600 x 4 - 83980 x 6 + 162945 x 8 - 192840 x 10 + 144876 x 12 
-70416 x 14 + 22170 x 16 - 4440 x 18 + 540 x 20 - 36 x 22 + x 24 ) 2 (4.41) 
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It is interesting to note that the second factor of Pjf 6 (x), (|4.41|) , is the same as -PfJ(x)/ 



' x 3 



which is the same polynomial appearing in ( 4.40Q and (|4.47|) . Again it should be stressed 



that these polynomials are monic and all the coefficients are integers. 

4.3.2 K = 56 for E 7 

Polynomial for 56, 

is slightly simplified for a different normalisation of /i: 

P^ 56 (x) = 2 28 Pj 7 6 (x/v / 2) = J] (x - /* • g) (A = V^/i, A 2 = 3) 

= 2044117922661550386613265625- 48583441852490416903125286500 x 2 

+ 403943437764362721049483097250 x 4 - 1594876299784237542505579618500 a; 6 

+ 3423181532874686547792360316875 x 8 - 4470973846715160163197028791000 x 10 

+ 3844463042762881314328636794900a; 12 - 2298706753677324429083230164600a; 14 

+ 994190889968661674517540390225 x 16 - 320292296385170629680242995500 x 18 

+ 78600569652362205629789205150 x 20 - 14948636823173617875192068460 x 22 

+ 2232949785098933644991402715 x 24 - 264680665744227895592493840 x 26 

+ 25089285771398909108223000 a; 28 - 1912398423761929885120080 a; 30 

+ 117632735062147883037411 x 32 - 5848529412061451267964 x 34 

+ 234966118304680273854 x 36 - 7609794291104570460 x 38 

+ 197734877929087065 x 40 - 4090765650038424 x 42 + 66612822142356 x 44 

- 839599815096 x 46 + 7991799795 a; 48 - 55327860 x 50 

+ 261954 x 52 - 756 a; 54 + a; 56 (4.43) 

4.3.3 K = A for E 7 and E 8 

The polynomials -P^ 7 (x) and -Pj^(x) are too long to be displayed here. See || for explicit 



1 /3 27 y<- 

forms. It should be stressed that five monic polynomials in x, P^f (x) (|4.40|) (and Pf J 



527, 



X 



(ggg)), P# 6 (x) ( ggg ), P^ 56 (x) ( |Q3D , P# 7 (x) and P# 8 (x) have integer coefficients only. 
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4.4 F 4 



The theory has two coupling constants gi and gs for the long (p 2 L = 2) and short (p| = 1) 
roots. We present the polynomials as a function of k = gs/ 9l- 

AAA 11 = A L for F A 

P 4 L (k\y) = P^(k\x) =\{{x-p.q)= J] {y-(p. q) 2 ) (p 2 L = 2) 

peA L peA L+ 

= 746496(1 + k) 6 (2 + k) 2 (l + 2k) - 4478976(1 + kf{2 + kf{l + 2k)y 
+124416(1 + kf{2 + kf{\ + 2ife)(91 + 64 k)y 2 
-13824(1 + k) 5 (2 + k) (2282 + 6049A; + 3712 k 2 + 512 A; 3 ) y 3 
+15552(1 + k) 4 (2 + k) (1718 + 5027A; + 4288 k 2 + 1024 A; 3 ) y 4 
-20736(1 + k)\2 + k) (695 + 1472A; + 704 k 2 ) y 5 
+1728(1 + A;) 3 (5878 + 16235A; + 14408A; 2 + 4096 A; 3 ) y 6 
-62208(1 + A;) 3 (38 + 71k + 32k 2 )y 7 + 432(1 + A;) 2 (838 + 1627A; + 784k 2 )y 8 
-576(1 + A;) 2 (62 + 61% 9 + 2160(1 + k) 2 y 10 - 72(1 + k)y u + y 12 . (4.44) 

4.4.2 11 = A s for F A 

P s A {k\y) EE P^(k\x) = J] (x - p ■ q) = J] (y-(p- q) 2 ) (p 2 = 1) 

peA s peA s + 
= 729A: 3 (1 + A;) 6 (2 + fc)(l + 2A;) 2 /4 - 2187A: 2 (1 + A;) 6 (2 + fc)(l + 2k) 2 y 

+243A;(1 + A;) 5 (2 + fc)(l + 2A;) 2 (64 + 9lk)/2y 2 

-27(1 + A;) 5 (l + 2k) (512 + 3712A; + 6049A; 2 + 2282A; 3 ) y 3 

+243(1 + A;) 4 (l + 2A;)(1024 + 4288A; + 5027A; 2 + 1718k 3 )/4y 4 

-162(1 + A;) 4 (l + 2k) (704 + 1472A; + 695A; 2 ) y 5 

+27(1 + A;) 3 (4096 + 14408A; + 16235A; 2 + 5878A; 3 ) y 6 

-1944(1 + A;) 3 (32 + 71k + 38A; 2 ) y 7 + 27(1 + A;) 2 (784 + 1627A: + 838A; 2 ) y 8 
-72(1 + A;) 2 (61 + Q2k)y 9 + 540(1 + k) 2 y 10 - 36(1 + k)y n + y 12 . (4.45) 

They are related with each other reflecting the self-duality of the F4 root system. If one 
replaces k by 1/k and y by y/{2k) in Pf(k\y), one obtains P[(k\y)/(2k) 12 : 

Pf(k\y) = (2k) 12 pS(l/k\y/2k), or P£ 4 L (k\x) = (2k) 12 P* s (1 / k\x / V2k) . (4.46) 
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It is well-known that F4 with the coupling ratio k = gs/gL = 2 is obtained from Eq by 
folding. This relates _F 4 polynomials to E$ polynomials: 

Pfr{2\x) = Pll{x)/x\ P£(2\x) {P^(2\x)) 2 = P£ 6 (x). (4.47) 



Both of them have trigonometric counterparts as will be shown later (|5.64 )- (p.66|) . The two 
polynomials Pp 4 L (k\x) and P^ As (k\x) have integer coefficients only. This property seems 
to be inherited from E e , too. 

4.5 G 2 

The theory has two coupling constants gi and g$ for the long (p| = 2) and short (p| = 2/3) 
roots. We present the polynomials as a function of k = gs/gL- 

4.5.1 U = A L for G 2 

P 2 L (k\y) = P£(k\x) =H(x-p.q)= J] (y~(p- if) (Pi = 2) 

peA L peA L+ 

= -27(1 + kf/2 + 81(1 + kf/A y - 9(1 + k)y 2 + y 3 . (4.48) 

4.5.2 K = A s for G 2 

P 2 s (k\y) = P A °(k\x) = J] (x - p ■ q) = J] (y-(p- qf) (p 2 = 2/3) 

= -k(l + k) 2 /2 + 9(1 + k) 2 /Ay -3(1 + k)y 2 + y 3 . (4.49) 

They are related with each other reflecting the self-duality of the G 2 root system: 

P 2 L (k\y) = (3k) 3 P 2 s (l/k\y/3k), or P^{k\x) = (3k) 3 P As (1 / k\x / V3k) . (4.50) 

The G 2 Calogero system with the coupling ratio k — gs/gL — 3 is obtained from that of _D 4 
by the three- fold folding Z) 4 — > G 2 . This implies analogous relations to ( 4.47|) 

P£(3\x) = Pl(x)/x 2 (K = V,S,S), P^(3\x){P^(3\x)) 3 = P A 4 (x). (4.51) 

Both of them have trigonometric counterparts, too, as will be shown later. The two polyno- 
mials Pq^ Al (k\x) and P^ As (k\x) have integer coefficients only. This property seems to be 
inherited from D4. 

Thirdly let us discuss the systems based on non-crystallographic root systems. 
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4.6 J 2 (m) 

The equilibrium points are easily obtained when parametrised by the two-dimensional polar 
coordinates |J: 

Q = (Qu = ^(sin <p, cos <p) , (4.52) 
-2 m 9 - TT _ 2 m(# e + # ) 

r = , 9? = (m : odd) ; r = , tan — \ — \jn '■ even), (4.53) 

a; 2m 2uo 2 \J g a 

in which g is the coupling constant in the simply laced odd m theory, whereas g (g e ) is the 

coupling constant for odd (even) roots in the non-simply laced even m theory. As 1Z we 

choose the set of the vertices of the regular m-gon R m on which the dihedral group him) 

acts: 

R m = |(cos(2j7r/m + to),sin(2j7r/m + t )) € M 2 j = 1, . . . , mj, (4.54) 
to = vr/2m (m : odd) ; t = (m : even). 

The polynomial n^ei?™ ( x ~ ^ ' 9*) * s obtained trivially: 

^(^) - = n & - ^ • 9) = n (* - s - (^f + ) ' (4 - 55) 

in which 

(po = ii (m : odd); ipo = 2 arctan Vk, k = g e /g (m : even). (4.56) 

For odd m P m (x) is proportional to the Chebyshev polynomial of the first kind T m (x) 
(see (|5.4|) ). For even m and for the equal coupling g e = g , P m (x) is also proportional 
to the Chebyshev polynomial T m (x) and thus the entire {P m (x) = 2 1 ~ m T m (x)} constitute 
orthogonal polynomials Jl|. For generic coupling g e 7^ g the orthogonality no longer holds. 
This can be seen most easily by the explicit forms of the lower members of P even in the 
non-singular limiting cases, g e = and g = 0: 



g e = : x 2 , x 2 (x 2 -l), x 2 (x 2 - 3/4) 2 , x 2 (x 2 - l/2) 2 (x 2 - 1), . . . , 

g o = : x 2 -l, (x 2 -l/2) 2 , (x 2 - l)(x 2 - 1/4) 2 , (x 4 - x 2 + 1/8) 2 , . . . , (4.57) 

which have definite sign in — 1 < x < 1. 

The following equivalences are well-known: A 2 = -^(3), B 2 = -^(4) and G2 = ^2(6). The 
^2(3) polynomial corresponds to the A 2 polynomial of vector V, 

P»foix) = ~T 3 (x) = ^=H 3 (V2x) = Pj^(x). (4.58) 

19 



As for ^2(4) polynomial, we obtain from (|4.55|) 



P% 4) (x)=x*-x 2 + k + , k = g e / 9o . (4.59) 
For the B 2 system, the Laguerre polynomial with a = k — 1 = g e j g — \ reads 

L ( 2 a \y) = \y 2 — (k+ l)y + k(l + k)/2, a = k-l. 
They are proportional to each other upon identification y = 2(l + k)x 2 . The 12(6) polynomial 



obtained from ( 4.55 ) reads, after some calculation 



Pf* (x) = x 6 - -x 4 + —x 2 t^—tt, k = gjg , (4.60) 

/ 2 (6)V ) 2 lg 16(1 + k) 1 y ' 



which is proportional to P 2 5 (^l?/) (H-49|) upon the same identification as above y = 2(1 + k)x 2 . 
4.7 H 3 and H A 

The non-crystallographic H3 and are simply laced root systems. In both cases the roots 
are normalised to 2, as with the other simply laced root systems, p 2 = 2. Then both monic 
polynomials Ph 3 ( x ) an d PhS x ) nave integer coefficients only. 

4.7.1 K = A for H 3 

P 3 A (V) = Ph 3 (x) = II (? ~ P ■ 9) = II (y ~ (P ■ 9) 2 ) (P 2 = 2 ) 
peA peA + 

= (-450 + 225y-30y 2 + ?/ 3 ) 

x (5625 - 22500 y + 27000 y 2 - 9600 y 3 + 1200 y 4 -60y 5 + y 6 ) 

x (22500 - 67500 y + 46 125 y 2 - 1 1700 y 3 + 1275 y 4 - 60 + y 6 ) . (4.61) 

4.7.2 ft = A for H 4 

p?(v) = Ph&) = U( x ~p-^= U (y-(p- 9) 2 ) (p 2 = 2 ) 

peA peA+ 

= (656100000000 - 1093500000000 y + 601425000000 y 2 - 154305000000 y 3 + 21343500000 y 4 
-1701000000 y 5 + 80392500 y 6 - 2250000 y 7 + 36000 y 8 - 300 y 9 + y 10 ) 

x (747338906250000000000 - 9964518750000000000000 y + 45172485000000000000000 y 2 
-90926233593750000000000 y 3 + 92928548278125000000000 y 4 
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-52841916742500000000000 y 5 + 18358385767875000000000 y 6 

-4169745135000000000000 y 7 + 648844128590625000000 y 8 - 71483472810000000000 y 9 
+5707114499700000000 y 10 - 335580296625000000 y 11 + 14683267406250000 y 12 
-480384270000000 y 13 + 1 1739694500000 y 14 - 212600700000 y 15 + 2804085000 y ie 
-26100000 y 17 + 162000 y 18 - 600 y 19 + y 20 ) 
x (1362025156640625000000000000000000 - 20430377349609375000000000000000000 y 
+ 110664543977050781250000000000000000?/ 2 - 280672524028933593750000000000000000?/ 3 
+406550434997274609375000000000000000 y 4 - 377089903500479578125000000000000000 y 5 
+240385775914964970703125000000000000 y 6 - 110467977515351929687500000000000000 y 7 
+37879740110299305937500000000000000 y 8 - 9947615045119592062500000000000000 y 9 
+2041289604408317542031250000000000 y 10 - 332506194678726581250000000000000 y 11 
+43529340095868749062500000000000 y 12 - 4624592400554729343750000000000 y 13 
+401746375286214215625000000000 y 14 - 28701181376029878750000000000 y 15 
+ 1693173350921514750000000000 y 16 - 82699244991680625000000000 y 17 
+3348318244893890625000000 y 18 - 112349936407545000000000 y 19 
+3118565868993450000000 y 20 - 71352951283125000000 y 21 + 1337980766062500000 y 22 
-20388872475000000 y 23 + 249452622000000 y 24 - 2408494500000 y 25 + 17897422500 y 2e 
-98550000 y 27 + 378000 y 28 - 900 y 29 + y 30 ). (4.62) 

5 Sutherland Systems 

Let us first discuss the systems based on the classical root systems. 

The equilibrium position is " equally-spaced" (see eq.(5.14) of [|l|) and translational invariant. 

r+l 

We choose the constant shift such that the coordinate of "center of mass" vanishes, qj = 0: 

j=i 

_ n(r + l-j) ixr _ ir ir(2j - 1) _ 
q >-^Ti WTT) ~ 2 " 2(r + l) " qr+2 -> V- 1 >-> r + 1 )- i 5 - 1 ) 
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5.1.1 TZ = V for A r 

For the vector weight fij G V ( f4.5|) , fij ■ q is independent on the constant shift of q. The 
above choice (15. If) leads to 



7T 7r(2j — 1) 7T 7T , . lN _ 

^■•g = 2~ 2(r+l) = % ~2 </V9< 2 U = 1 '--- r + 1 )- ( 5 - 2 ) 



In this case the polynomial (O) is given by 

r+1 r+1 



P r v (x) = PX» 



f[ (x - frinfo ■ §)) = f[ fx - cos ffi 1} ) = 2-T r+1 (x) . (5.3) 
i=i i=i V I + J / 

Here T n (cos<^) = cos(n^) is the Chebyshev polynomial of the first kind, whose Rodrigues' 
formula is 

T n ( X ) = - ^^". (l - X 2 ) 1 / 2 f AV(1_ x 2)n-l/2 = 2 n-l x n + . . . (54) 



(2n-l)!! v 
They are orthogonal to each other: 



1 P r v (x)P s v (x) , . . . 

7== oc £ rs . (5.5) 
-i V 1 — £ 

This is a new result. Another definition 

r+1 

= - 2sin(^- • q)) = 2T r+1 (x/2) = 2 r+1 P^ s (x/2) (5.6) 

provides a monic polynomial with a// integer coefficients. 
It is easy to see that 

^X.cfa) = Il(z - cosC/i, • g)) = Y[ (x - sin ^ + 1) ) 

does not give rational polynomials, for example, PY Jx) = x 2 — V^x + 1/2. In fact, in most 
cases the polynomial Pj c {x) is not of rational coefficients. In the rest of this paper we will 
not consider this type of polynomials. 
The other polynomials, 

r+l r+l 



pla*) = n(* - sin ( 2 ^- • *)) = n (* - sin 

i=i ' j=i V ' 

p A r ,c2^) = n ( x ~ c ° s ( 2 / i i • ?)) = n ( x + c ° s 71 /, x ) 

.7=1 .7=1 ^ ' 



J=l 
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are essentially the same as P][ r s (x), ( |5.3|) . Only the constant term can be different: 



777" 

PlA x )- p lM = -2- r sin- P% tta {x) -Pl,.(x) = (-1)^2- (5.7) 



Thus we consider only the polynomial P^ r s (x) = Y\ ( x ~ sin(/i • q)\ for various 1Z of A r . 



5.1.2 K = Vj for A r 

From (|4.9|) and ( |5.2|) , the polynomial ( p.4|) is given by 

P X» = II { x - sin fe + ■ • '&)) = P X^ 1 " i (^)- (5-8) 

l<ji<-<j<<r+l 

This polynomial can be expressed in terms of the coefficients of T r+1 (x) by the same method 
as given m section gX§ and 2 Di Pj; s (x/2) gives a monic polynomial with integer coefficients. 
See || for the explicit forms of the polynomials P^* s (x) of lower rank r. 
As in the Calogero case, we report only on V2: 



>v 2 



Y[ \x- sin(g j + 

l<j<i<r+l 



(5.9) 
r : odd 



l<j<Z<(r+l)/2 



X 



r/2 



r/2 

J^J(x 2 — sin 2 gj) ■ jQ ^x 2 — sin 2 (g.,- — q~i)j (x 2 — sin 2 (gj + q~i) J r : even. 



3=1 i<i<z<r/2 

Based on the fact that the zeros of Chebyshev and Jacobi polynomials are related as seen from 
the formulae (|5.28|) and (|5.29|) , this can be expressed by using the polynomials associated 
with the BC r Sutherland systems in the following way: 



PaLiM = ^ rir ~ l) ^ P BC r U^ 1/2|1 - 2x 2 



P 



V 2 1 



2-^x r - 1 Pl r Jx)P^ c2 (l, 1/2|1 - 2x 2 ). 



(5.10) 
(5.11) 



The explicit forms of the functions Pbc^c2{^i for lower r are given in section p. 2. 5 



5.1.3 K = A for A r 

The polynomial has a factorized form: 



IJ (s 2 -sin 2 (g j -g i ; 



l<j<Z<r+l 



x — 1 



(r = l) 
(r > 2) 



(5.12) 



23 



It is elementary to evaluate P£ s (x) for lower rank: 

pim= n (« , --» , (^)) 

l<j<Z<r+l V 7 

= 2" 6 (4x 2 -3) 3 , 

PtM = 2~ V - l) 2 (2x 2 - l) 4 , 

P^ )S (x) = 2- 20 (5 - 20x 2 + 16x 4 ) 5 , 

Pi g ( x ) = 2 - 2 V - l) 3 (4x 2 - l) 6 (4x 2 - 3) c 



P Ae,s( x ) = 2~ 42 (-7 + 56x 2 - 112x 4 + 64x 6 ) 7 . 

For r = 1, 3 and 5, Pj£ s (x) are of definite sign in — 1 < x < 1. They can never be orthogonal 
with each other for whichever choice of the positive definite weight function. 

5.2 BC r and D r 

As shown in |]], the equations ( |2.7|) for A = £>C r read 

^> sin2g 7 - cosg 7 - n cos2g 7 - ^ 10 \ 

-2(7m> ■ + ^ + 2^ L — = j = l,...,r. 5.13 

cos 2g, — cos 2qi sin sm 2^- 

For non- vanishing and gi, sin 2gj = cannot satisfy the above equation. Thus by dividing 
by sin 2^ we obtain for k x = gs/gu, h = gi/gM- 



1=1 X 3 



in which Xj = cos2qj. These are the equations satisfied by the zeros {xj \j = 1, . . . ,r} of 
the Jacobi polynomial Pr a,l3 \x) f| with 

a = k 1 + k 2 -l, (3 = k 2 -l. (5.15) 

The Rodrigues' formula for the Jacobi polynomial P^" ,l3 \x) reads 

P^\x) = t^(l-x)~ Q (l + x)-^ (-f) ( (1 - x) n+a (l + x) n+ ?) (5.16) 



2 n n\ \dx 
1 T(2n + a + (3 + l) n 
2 n n\ T(n + a + p + 1) X + 
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For A = D r , we have gs = 9l = 0, implying a = (3 = — 1. We choose 
qi = 0, q r = 7r/2 ( < ^=>- cos2gx = 1, cos2g r = —1), 

then fl2~7f ) read 

E^- + -^T + ^T = (J=2,...,r-1), (5.17) 

in which 2j = cos2qj (j = 2, ...,r — 1). These are the equations satisfied by the zeros 
{xj | j = 2, . . . , r — 1} of the Jacobi polynomial Pj^! 2 (x) ||. In fact, there is an identity 

4P r (- 1 '" 1 )(x) = (x 2 - l)P r ( S(x), (5.18) 

which means that {l,x 2 , . . . , aV-i? —1} are the zeros of Pr~ X ~ (x). This allows to treat D r 
as a limiting case of BC r . 

The possible 7?.'s for BC r are A5, A M and A L . Since A5 = {±ej \j = 1, . . . ,r} and 
A^ = {±2ej \ j — 1, . . . , r}, we have trivial identities among the polynomials 

Pncrjki, h\x) = P^S r , s2 (h, k 2 \x), P££ r>c (h, k 2 \x) = P*3 r ,c2(ku k 2 \x). (5.19) 

In other words, these relations prompted us to introduce the polynomials of the forms 
Il/iS-ft ( x — si n (2/i • q)) and n^erc ( x — cos(2// • qj). For BC r Sutherland system we consider 
1Z = As and A m only. 

5.2.1 K = A s for BC r 

Since As is even and that {xj = cos2qj \j = 1, . . . , r} are the zeros of the Jacobi polynomial, 
it is natural to consider the polynomial (|3.8| ) 

Pb^A^ « = II (* - cos2 ^ = 2V! r(2r+l+^ + l) Pr(Q,/3)(x) ' (5 - 20) 
with a = fci + A; 2 — 1 and (3 = k 2 — 1. They are orthogonal to each other: 

j P^\x)P^\x)(l-x) a (l + xY dx oc5 rs . (5.21) 

As remarked in (|3l9|) , the polynomial Pbc + &{ki, k 2 \x) is equivalent with P^? s (fc 1? ^l^)- 
Needless to say that 2 n n\Pn°'' l3 \x) is a polynomial in the parameters a and /3 with integer 
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coefficients. Thus P B q + c2 (ki, k 2 \x) ( |5.20|) is a rational function in a and (3 with integer 
coefficients. 

The other polynomial P B q s2 {ki, k 2 \x) can be easily obtained by (|3.10|) : 



r 

3=1 



in which u 2 = 1 — x 2 . Remark that P r 



5.2.2 K = V for D r 



-X) 



-iyp r w ' a \x). 



This is a special (ki = k 2 = or a = (3 = — 1) case of the previous example. As in the 
previous example, we introduce 

P?+( X ) ee Pj+ c2 (x) = J! (x - cos2g,) = 2 -^^Py i '- 1 \x) 



(x 



3=1 
r-l 

+i)( x -i)H( x - 

3=2 



' (2r-2)! 



(2r - 2)! 



-(x+l)(x- l)P^ 2 '(x).(5.23) 



They are orthogonal to each other: 



P r v+ 0) P7 + (x) ( 1 - x) - 1 ( 1 + x) - 1 dx oc y P r ( i'J } (z) P^? (x) (1 - a:) ( 1 + x) dx oc5 rs . 

(5.24) 

Corresponding to the Dynkin diagram folding D T+ \ — > B r and ( |5.23 ), we obtain 



(* - l)^£ c2 (2, 0|x) = P^+ +1 , c2 (*) = P^ + +ljC2 (0,0|x) 



A S 4 



(5.25) 



which is the trigonometric counterpart of ( |4.21| ). 

The other polynomial Pp s2 ( ;r ) nas a simple form 



r 

= n( x2 - sin2 ( 2 ^) 

3=1 



■1) 



P ^V!(r-2)r, 



(5.26) 



u 2 — >1— a; 2 



(2r - 2)! 

which is of a definite sign in — 1 < x < 1. Thus they do not form any orthogonal polynomials. 
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5.2.3 A-ir-i C r and the relationship between Chebyshev and Jacobi polyno- 
mials 

As in the Calogero case, the Dynkin diagram folding A2 r ~i — > C r implies 

PL-U*) = (-2)- r ^: c2 (0, 1/2|1 - 2x 2 ). (5.27) 
Indeed there are relations between Chebyshev and Jacobi polynomials: 

2^T 2r (x) = (-l)^ ^i| p(-V2,-i/2) (1 _ 2x 2 ); (5>2g) 

2~ 2 'T 2r+1 (x) = (-l)''g|xP( 1 / 2 '- 1 / 2 )(l-2x 2 ), (5.29) 

on top of the well-known relations (eq(4.1.7) of ||): 

1 ■ 3 ■ ■ • (2r - 1 



T r {x)=Py i ' 2 ^ 2 \x). 



2 • 4 • • • 2r 

The former corresponds to (5.27) and the latter implies 



Pj 2r » = (-2)-a;P^ c2 (l, 1/2|1 - 2x 2 ). (5.30) 
5.2.4 ft = S and S for D r 



As in the Calogero systems, the symmetry of the D r Dynkin diagram implies that Pg a (x) = 
P§ r a (x), a = s, c, s2, c2. Among them P^, S c (x) do not always give rational polynomials. As 
remarked above ( |3.9|) , P^ s s (x) are equivalent to Pnt'c2 + ( x ) f° r even rank r. Thus we list for 
lower rank r the polynomials Pr^'c2 + ( x ) an< ^ ^n^si (• r ) : 



>V+,S + ,S+/ \ 
D 4 ,c2 W 


= (x 2 - 


l)(x 2 -l/5), 




(5.31) 




4/ 2 

= x (x 


-4/5) 2 , 




(5.32) 




= (x 2 - 


l/2) 4 (x 4 -x 2 - 


1/196) 2 , 


(5.33) 


-^D5,s2( X ) 


= (x 2 - 


l/2) 4 (x 4 -x 2 - 


1/196) 2 , 


(5.34) 


pS+, S-j- / \ 
- r D 6 ,c2 W 


= 3~ 4 7" 


3 x 4 (21x 4 - 28x 2 


+ 8) 2 (63x 4 - 72x 2 + 16), 


(5.35) 




= 3- 8 7- 


6 (x 2 - l) 4 (21x 4 


- 14x 2 + l) 4 (63x 4 - 54x 2 + 7) 2 . 


(5.36) 



It is interesting to note that the formula ( p.lOD applies to D 5 (conjugate) spinor represen- 
tation S (S), which is not even. This is because the set of values {fi ■ q \ fi G S} is even. 
Moreover, the function in Q5.33|) is invariant under x 2 — > 1 — x 2 . 
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5.2.5 K = A M for BC r 

The set of middle roots is A M = {±(e,- — e/),±(ej + e/) | 1 < j < I < r}. As in the 
Calogero systems in |4.2.5| , the polynomial Pg™ s (k\x) can be expressed neatly in terms of the 
coefficients of the polynomial Pg S s (k\x). Suppose we have two polynomials in y: 



f = f[(y - sin 2 x t ) = £(-l) 4 a, V n ~\ (5-37) 

i=l i=0 

9= \\ (y - sin 2 (xi- Xj))(y - sin 2 ^ + Xj)\ (5.38) 



l<i<j<n 

Let us denote bi = sin 2 a;,, then we obtain g as a symmetric polynomial in bf 



9= 11 (y 2 -2(b l + b J -2b i b J )y + (h-b J ) 2 ^ (5.39) 

l<i<j<n 

and {ai} are the basis of the symmetric polynomials in bf. 

a i= Yl h h---K- ( 5 - 4 °) 

l<ii<— <ji<n 

Thus g can be expressed in terms of the coefficients {a{\ of / with integer coefficients. For 
example: 

n = 2 : g = y 2 - 2(a x - 2a 2 )y + a\ - 4a 2 , (5.41) 
n = 3 : 

g = y 6 — 4(ai — a 2 )y 5 + 2(3a 2 — a 2 — Aa\a 2 — 12a3 + 8aias)y 4 

— 2{2a\ — a\a 2 — 13a 3 — 2a 2 a 2 — 4a| — 2aia 3 + 32a 2 a 3 — 32a 2 )?/ 3 

, a o o 0\ O 

+ (a x + 2(2^2 — 7a 2 — 2Aa\a 3 — 8aia 2 — 16a 1 a3 + 120a2£i3 + 16aia2a3 — 144a 3 )y 

,o 00 0\ 

— 2{a\a 2 — 3aia 2 — 90^3 + 27a 2 a 3 — 2a 2 — 2a^a 3 + 18aia 2 a 3 — 54a 3 )y 

+ a\a\ - 4a 2 - 4a 3 a 3 + 18aia 2 o.3 - 27a 2 . (5.42) 

If / is of rational coefficients, so is g. 

Here are some explicit forms of P b q t s (ki, k 2 \x) for lower rank r (see also ||): 

p a m ( . Ux) 4(l + k 2 )(l + k 1 + k 2 ) 4(l + k 2 )(l + k 1 + k 2 ) 

21 j " (1 + *i + 2fc 2 )(2 + h + 2k 2 f (1 + h + 2^ 2 )(2 + ft + 2fc 2 ) y+y ' l5 ' 4dj 
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PBC 3 ,s(h,h\x) 



108(1 + fc 2 )(2 + k 2 ) 2 (l + h + k 2 ){2 + fci + k 2 ) 2 
(2 + h + 2k 2 ) 2 (3 + £4 + 2fc 2 ) 3 (4 + fci + 2A; 2 ) 4 
108(1 + fc 2 )(2 + fc 2 ) 2 (l + fci + fc 2 )(2 + fei + fc 2 ) 2 (10 + 3fci + 6fc 2 ) 
(2 + h + 2fc 2 ) 2 (3 + h + 2A; 2 ) 3 (4 + ^ + 2A; 2 ) 4 



9(2 + £; 2 ) 2 (2 + k t + k 2 ) 



(164 + 196/si + 41fc 2 



(2 + fci + 2A; 2 ) 2 (3 + ki + 2k 2 ) 2 (A + h + 2k 2 y 

+392A; 2 + 292/^^2 + 32fc 2 A; 2 + 292A; 2 + 96hkj + 64Jbf) y 2 

— ; 4( i +fc2 l (2 + " 1 ^" 2) i r^s (792 + 1278^ + 639fc 2 

(2 + ki + 2k 2 ) 2 (3 + h + 2k 2 ) 2 (A + h + 2k 2 ) 3 K 1 

+99Jfe? + 2556A; 2 + 3088^^2 + 1052/c 2 /c 2 + 88kfk 2 + 3088fc| 

+2562ife a A;| + 504fc 2 fc 2 + 16fc 3 fc 2 + 1708^ + 832A; X ^ + 64A; 2 fcj 

+416A; 2 + 80fcifc^ + 32A|) y 3 
6(2 + fc 2 )(2 + kx + A; 2 )(26 + 17Jfei + 34A; 2 + 8^^ + 8A; 2 
(2 + h + 2k 2 )(3 + h + 2k 2 )(A + h + 2k 2 ) 2 
12(2 + A; 2 )(2 + A; 1 + A; 2 ) 5 6 

-y + y ■ 



-y 



(3 + /c 1 + 2/e 2 )(4 + A; 1 + 2A; 2 ; 
5.2.6 TZ = A for A 

These are the fci = and k 2 = limit of the formulae given in the previous subsection 



(5.44) 



pA+ 




= (l + x) 3 (-3/5 + x)(-l/5 + x 2 ) 4 , 


(5.45) 






= a; 6 (-4/5 + a; 2 ) 8 (-16/25 + x 2 ), 


(5.46) 


pA+ 




= x 4 (l + x) 3 (-l/7 + x)(-4/7 + x 2 ) 2 (-3/7 + x 2 ) 4 , 


(5.47) 


Pd 5 ,s2 


» 


= (-1 + x 2 ) 4 x 6 (-4/7 + x 2 ) 8 (-3/7 + x 2 ) 4 (-48/49 + x 2 ), 


(5.48) 


pA+ 

- r D 6 ,c2 




= 3~ 9 7~ 7 (1 + x) 4 (-3 - 14x + 21x 2 )(l - 14x 2 + 21x 4 ) 4 (7 


- 54x 2 + 63x 4 ) 2 , (5.49) 


Pd 6 ,s2 




= 3~ 18 7~ 14 x 8 (8 - 28x 2 + 21x 4 ) 8 (16 - 72x 2 + 63x 4 ) 4 (128 - 


- 560x 2 + 441x 4 ).(5.50) 



It is trivial to verify that ( |3.10| ) are satisfied: 
The Dynkin diagram folding D T+ i — > B r relates the functions 



-r c 



BCr) j2,0\ X ) (P£;j2,0\ X )) = P£ +1 Jx) = 
which is the trigonometric counterpart of the identity ( [4.37 ). 



^: + +1 , c2 (o,oix), 



(5.51) 



(5.52) 



Next we discuss the systems based on the exceptional root systems. As in the Calogero 
systems, we have relied on numerical evaluation of the equilibrium points. 
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5.3 E r 

5.3.1 K = 27 and A for E 6 

We have evaluated two polynomials independently: 
P eIc2(%) = II ( x ~ cos(2/i • q)^ 



[-l +x) 3 (l +2xf 

2 18 7 4 11 6 



(-743 - 42651 x + 708939 x 2 - 1704045 x 3 - 1890504 x 4 



+7043652 x 5 + 1260336 x 6 - 9391536 x 7 + 4174016 x 9 ) , (5.53) 

and 

P eIs2( x ) = Yl (x - sin(2fi ■ q) ^ 



x 3 (-3 + 4x 2 ) 3 



(-221709312 + 39409774992 x 2 - 786312492840 x 



4 



2 18 7 4 U 6 

+6804048466593 x 6 - 32072860850184 x 8 + 89147361696624 x 10 
-149154571577088 x 12 + 147001580732160 x 14 - 78400843057152 x 16 
+ 17422409568256 x 18 ) . (5.54) 

Although the set of minimal weights 27 is not even, that is —27 = 27 ^ 27, these two 
polynomials are related. The "formula ( 3.10]) " is valid, 



-u 



(5.55) 

i 2 ^l — X 2 



This is the same situation encountered in the D 5 (conjugate) spinor representations S (S) 
in ( |5.33| ). This provides a strong support for the above results. 
As for 1Z = A, we have: 



peA H 

(2x+ l) 6 , , o N 

2 V 24 ; n (-235 - 627 x + 231 x 2 + 847 x 3 ) 

x (-743 - 42651 x + 708939 x 2 - 1704045 x 3 - 1890504 x 4 

+7043652 x 5 + 1260336 x 6 - 9391536 x 7 + 4174016 x 9 ) 3 , (5.56) 



Pk,J?) = II (^-sm(2p-g) 



peA H 



30 



f_3 + Ax 2 ) 6 

y Ta ^ — (-48384 + 422928 x 2 - 1036728 x 4 + 717409 x 6 ) (5.57) 
x (-221709312 + 39409774992 x 2 - 786312492840 x A + 6804048466593 x 6 
-32072860850184 x 8 + 89147361696624 x 10 - 149154571577088 x 12 
+147001580732160 x 14 - 78400843057152 x 16 + 17422409568256 x 18 ) 3 . 



5.3.2 TZ = 56 for E 7 

We have evaluated two polynomials independently: 
P ErM x "> = II (^-cos(2^-g) 



4 

X (9332954265600 - 345319307827200 x 2 + 5422446428313600 x 4 

J. J. -Lo .Li 

-47902580312348160 x 6 + 266584469614182720 x 8 
-991356255189780480 x 10 + 2543382104409514368 x 12 

16 



-4564307435286703104 x 14 + 5717674981551733200 x 
-4899020276961851040 x 18 + 2736363552042360240 x 20 
-897719270582318184 x 22 + 131214258464743597 x 24 ) , (5.5* 



and 



peUw = n { x - sin ^ • s)) 

6 

-1 +X 



„2N4 



(7824285157 - 1019921980260 x 2 + 44927774191218 x 4 



H8 13 10 17 12 

-933762748148260 x 6 + 10512912980210355 x 8 
-70729109671077000 x 10 + 302444017343367900 x 12 
-850322103495681960 x 14 + 1590230624766864795 x 16 



+957192223677842580 x 18 + 1521592634309937618 



x 



20 



-676851830994604980 x 22 + 131214258464743597 x 24 ) 2 . (5.59) 
These two polynomials satisfy ( |3.10| ), 



(5.60) 
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5.3.3 K = A for E 7 and E h 

,A 



The polynomials P^ s r 2 )( x )i Pe^c2( x )i r = 7, 8 are too long to be displayed here. Their 
degrees are 63 and 126 for E<j and 120 and 240 for E 8 . They are given in ||. They all satisfy 
the consistency condition (|3.10| ) 



Pl,M = P£MP£*{-») 2 2 (r = 6,7,8). (5.61) 
at the level of each factor. 

5.4 F 4 

We present the polynomials as a function of k = gs/ 9l- The polynomials Pp^ 2 ' As+ (k\x) and 
P^ i ^ s '(/c|x), satisfying the condition (|3.1Q) , are too lengthy to be displayed here. They are 
given in ||. Here we give Pp 4 L ' s As (k\x) which have shorter forms. As before we use y = x 2 . 

5.4.1 K = A L for F 4 

p£ s (k\y) = p^ s (k\x) = n { x - sin (p ■ ?)) = n (y - sin2 (p ■ «)) 

peA L peA i+ 

2 12 3 6 (1 + kf(2 + £Q 2 (3 + £Q 3 (1 + 2k) 

~ (3 + 2fc) 3 (4 + 3fc) 4 (5 + 3fc) 5 (6 + 5A:) 6 

2 13 3 6 (1 + kf{2 + k) 2 {3 + A;) 3 (l + 2k){U + 9k) 

(3 + 2£;) 3 (4 + 3A;) 4 (5 + 3k) 5 {6 + 5fc) 6 V 

2 U 3 6 (1 + fc) 5 (2 + fc) 2 (3 + A;) 3 (l + 2A;)(232 + 346A; + 123A; 2 ) 2 

+ (3 + 2fc) 2 (4 + 3fc) 4 (5 + 3/b) 5 (6 + 5A;) 6 ^ 

_ 2 n 3 4 (l + fc) 5 (2 + k)(3 + kf 

~~ (3 + 2A;) 2 (4 + 3A:) 4 (5 + 3fc) 5 (6 + hkf 

x (30432 + 133672A; + 211560A: 2 + 155726A; 3 + 54075A; 4 + 7128A; 5 ) y 3 

2 8 3 6 (1 + fe) 4 (2 + fc)(3 + fc) 2 

+ (3 + 2A;) 2 (4 + 3/fc) 3 (5 + 3fc) 4 (6 + 5fc) 6 

x (19296 + 90360fc + 159652& 2 + 137582A; 3 + 61155A; 4 + 13264fc 5 + 1088fc 6 ) y 4 

2 9 3 4 (1 + k) 4 (2 + k)(3 + k) 2 , , 

-7 r-T^ ttVt —7—7- — - — — (283824 + 1395972& + 2711556A; 2 

(3 + 2A;) 2 (4 + 3k) 3 (5 + 3k) 4 (6 + 5k) 6 v 

+2704381A; 3 + 1489217A; 4 + 447066A; 5 + 65952A; 6 + 3456A; 7 ) 

2 7 3 4 fl + k) 3 (3 + k) 2 

+ 7 1NO/ v 1NQ ; v -4t — (1046592 + 6283632A; + 15907184A; 2 

(3 + 2fc) 2 (4 + 3fc) 3 (5 + 3A;) 3 (6 + 5A;) 6 v 

+22205264A; 3 + 18708264A; 4 + 9754573A; 5 + 3088726A; 6 



y 5 
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+553392A; 7 + 47232A; 8 + 1152A; 9 ) y 
2 8 3 4 (l + k) 3 (3 + k) 2 



(35736 + 163412A; + 300546A; 2 



(3 + 2A;) 2 (4 + 3k) 2 (5 + 3fc) 3 (6 + 5A;) 5 

+286499A; 3 + 151260A; 4 + 43412A; 5 + 6048A; 6 + 288A; 7 ) y 7 
864(1 + A;) 2 (3 + k) 



(3 + 2k)(4 + 3/fc) 2 (5 + 3A;) 2 (6 + 5A;) 4 

x (33120 + 130392A; + 199564A; 2 + 150034A; 3 + 57649A; 4 + 10632A; 5 + 720A; 6 ) y & 
1152(1 + A;) 2 (3 + k) 



3 

9 



(3 + 2ife)(4 + 3fc) 2 (5 + 3A;) 2 (6 + 5k) 

x (3312 + 10668A; + 12946A; 2 + 7313A; 3 + 1899A; 4 + 180A; 5 ) y 
144(1 + A;) 2 (3 + A;)(116 + 133A; + 30A; 2 ) 10 
+ (3 + 2k) {4 + 3k) {5 + 3k) {6 + 5k) 2 V 
72{l + k){3 + k) n , 12 



(5 + 3A;)(6 + 5A;) 
5.4.2 11 = A s for Fa 



y Ll + y 12 . (5.62) 



P 4 S s (k\y) = P*%(k\x) = J] (x - sin(p • g)) = J] (y - sin 2 (p • q)) 

peA s peA s+ 

729A; 3 (1 + k) e (2 + k)(3 + k)(l + 2k) 2 
(3 + 2A;) 2 (4 + 3A;) 3 (5 + 3k) 3 {6 + 5kf 
2916A; 2 (1 + kf{2 + fc)(3 + fc)(l + 2fc) 2 (9 + 7 k) 

(3 + 2A;) 2 (4 + 3fc) 3 (5 + 3A;) 3 (6 + 5kf V 
1458Jfe(l + kf(2 + fc)(3 + fc)(l + 2A;) 2 (48 + 115k + 58A; 2 ) 2 
+ (3 + 2A;) 2 (4 + 3A;) 2 (5 + 3kf{Q + 5A;) 5 V 

324(1 + kf(3 + k)(l + 2k) 
~ (3 + 2A;) 2 (4 + 3k) 2 {5 + 3A;) 3 (6 + 5A;) 5 

x (1152 + 11712A; + 33125A; 2 + 38811A; 3 + 20104A; 4 + 3804A; 5 ) y 3 
729(1 + A;) 4 (l + 2A;)(1536 + 8960A; + 17519A; 2 + 15049A; 3 + 5788A; 4 + 804A; 5 ) 4 
+ (3 + 2jb)(4 + 3A;) 2 (5 + 3fc) 3 (6 + 5k) 4 V 

324(1 + fc) 4 (l + 2fc) 

~~ (3 + 2l)(4 + 3A;) 2 (5 + 3fc) 3 (6 + 5A;) 4 

x (26496 + 112704A; + 177478A; 2 + 130823A; 3 + 45354A; 4 + 5913A; 5 ) y 5 
162(1 + A;) 3 
+ (3 + 2k){4 + 3A;) 2 (5 + 3kf{6 + 5A;) 3 

x (37824 + 208304& + 455436A; 2 + 505691A; 3 + 300828A; 4 + 90935A; 5 + 10902A; 6 ) y 6 
324(1 + A;) 3 (9984 + 42832A; + 70360A; 2 + 55311A; 3 + 20783A; 4 + 2978A; 5 ) 7 
(3 + 2Ai)(4 + 3k) 2 {5 + 3A;) 2 (6 + 5A;) 3 y 



33 



54(1 + fc) 2 (4224 + 13765A; + 16027A; 2 + 7876A; 3 + 1380A; 4 ) 8 
+ (3 + 2ife)(4 + 3fc)(5 + 3fc) 2 (6 + 5k) 2 V 

72(1 + A;) 2 (2 + A;)(345 + 628A; + 276A; 2 ) 9 



(3 + 2Jfe)(4 + 3fc)(5 + 3fc)(6 + 5fc) J 



' 2/ 



36(l + fc) 2 (52 + 29fc) _ 36(1 + A;) 
(4 + 3fc)(5 + 3fc)(6 + 5fcr (6 + 5fc) y y 1 ; 

The folding P 6 — » P 4 relates E$ polynomials to P 4 polynomials at the coupling ratio k = 
9s/ 9l = 2. We have corresponding to ( |4.47| ) 

P£*(2\x) = Pll s2 (x)/x\ P^ c2 (2\x) = Pil c2 (x)/(x - l) 3 , (5.64) 
P£ a (2\x) (P£f a (2|x)) a = P*Jx) (a = s, s2), (5.65) 

(p£%mx)) 2 = p£a*)- ( 5 - 66 ) 

The self-duality of the P 4 Dynkin diagram relates A^ polynomials to As ones. For example, 
we obtain: 



mpf a (2\ y ) ^PL{2\y) 



(5.67) 
(5.68) 
(5.69) 



847y 3 - 1386y 2 + 594y - 27 {Ay - 3) 3 ' 

717409 P 4 L s2 (2| j/ ) _ 64P 4 g s2 (2| ?/ ) 

717409?/ 3 - 1036728?/ 2 + 422928?/ - 48384 {Ay - 3) 3 ' 

U7P^ 2 {2\x) 8P^+(2|x) 
847a; 3 + 231x 2 - 627a; - 235 {2x + l) 3 ' 

which are a factor of the parent polynomials Pj^, s , Pj^ s2 and Pe 6 + c2 , respectively. 

5.5 G 2 

Two types of polynomials n p erc + ( x ~ cos (2p • q)) and Ylpen ( x — sin(2p • q)) are evaluated. 
For the latter we use y = x 2 . 



5.5.1 U = Ar for G 2 



Pq?Mz) = J] (^-cos(2p-g) 



P6A i+ 

27 - 81fc - 99A; 2 + 107A; 3 + 80 k A - 16 k 5 3(27 - 81 k 2 - 40 k 3 + 16 A; 4 ) 
2(2 + k) 2 (3 + 2 A;) 3 + 2(2 + A;) (3 + 2 A;) 3 X 

3(3 + 2 fc -2 A; 2 ) „ , . . 
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PtMv) = Pc 2 L Ak\x)= n fx - sin(2p • g)) = f] (y - sin 2 (2p ■ q) 

peA L 

729 (1 + k) 2 (-3 + A; + 8 A; 2 ) 2 



4 (2 + A;) 4 (3 + 2 A;) 5 
729 (1 + A;) 2 (6 + 13 k + 8 A; 2 ) (9 - 6 k + 13 A; 2 + 8 A: 3 ) 

4(2 + A;) 3 (3 + 2 A;) 6 V 
27 (1 + A;) (9 + 12 A: + 13 A; 2 + 8 A; 3 ) , , 



(2 + A;) 2 (3 + 2 A;) 3 



(5.71) 



5.5.2 ^ = A S for Go 



^x — cos(2p ■ g) 



9 — 21 fc — 13 A; 2 + A; 3 3(-3 - 4 A; + A; 2 ) 3 A; „ „ .„ „„. 

ac + „ — ; x 2 + x 3 , (5.72) 



2(2 + A;)(3 + 2A:) 2 2(2 + A:)(3 + 2 fc) 3 + 2 fc 
= PgIsM*)= II - sm(2p • g)) = J] (y - sin 2 (2p ■ g)) 
27 (-3 + A;) 2 jfc (1 + A;) 2 27 (1 + A;) 2 (9 + 12 k + A; 2 + 2 A; 3 ) 



4 (2 + A:) (3 + 2 A;) 4 
9 (1 + k) (3 + 2 k + A; 2 ) „ , 



4(2 + A;) 2 (3 + 2A;) C 



(2 + ife) (3 + 2 fc) 2 
They satisfy the formula ( |3.10| ) 



(5.73) 



p, 



Aw 



G 2 ,c2 



G 2 ,c2V 



(.:2 V 



-u 



i 2 —*l—x 2 

(5.74) 



The Dynkin diagram folding P4 — > G 2 implies 

^(Slx) = Pfcjx)/{x - 1), P^ s2 (3|x) = P£ 4 , s2 (x)/x 2 (K = V,S,S), (5.75) 
^.(3|«) (PgU 3 \ x )T = P dU*) (a = *, «2), (5.76) 



^: 2 (3k) Pg 2 :: 2 (3| 



x. 



p, 



A+ 



D 4 ,c2 



(5.77) 



which correspond to ( (4.51| ). The self-duality of the G 2 Dynkin diagram relates P G £ s ( 2 )(3|x) 



and P£% (2) (3|x) (see § for P£f s {k ,•)): 



5P 2 L s (3|y) F*,(3|i/) 25P^ 2 (3|y) P 2 % 2 (3|y) 



5y-l 
5P 



g^(3|x) 
5x — 3 



y - 1 ' 25y - 16 

^(3|x) 



x + 1 
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(5.78) 
(5.79) 



which are a factor of the parent polynomials, Pq 4 a , P^ 4 S 2 an d ^ot c2-> respectively. 

6 Summary and Comments 

We have derived Coxeter (Weyl) invariant polynomials associated with equilibrium points 
in Calogero and Sutherland systems based on all root systems. For the classical root sys- 
tems, the polynomials are well-known classical orthogonal polynomials; Hermite, Laguerre, 
Chebyshev and Jacobi of the degree equal to the rank r of the root system (r + 1 for the 
A r case), when the smallest set of weights 1Z are chosen. For the other choices of 72-'s, 
the polynomials are related with the corresponding classical polynomials but they no longer 
form an orthogonal set. For the exceptional and non-crystallographic root systems, these 
polynomials are new. Some polynomials are given in |J, since they are too lengthy to be 
displayed in this paper. These new polynomials have (much) higher degree than the rank r; 
27 and 36 for E 6 , 28 and 63 for E 7 , 120 for E 8 , 12 for F 4 , 3 for G 2 , m for J 2 (m), 15 for H 3 
and 60 for if 4 . Defined only for sporadic degrees, these new polynomials do not have the 
orthogonality property, except for those corresponding to the dihedral group l2{m) with the 
uniform coupling g = g e = g Q . In this case Chebyshev polynomials are obtained [|TJ. 

All these new polynomials share one remarkable property with the classical polynomials; 
their coefficients are rational functions of the ratio of the coupling constants with all integer 
coefficients. In most cases, they are monic polynomials with integer coefficients only. 

In the rest of this section, we give a heuristic argument for "deriving" the classical 
orthogonal polynomials with the proper weight function from the pre-potential W ( |2.4| ) at 
equilibrium. We add one degree of freedom, a new coordinate q r+ i (q r +2 for A r ), to the rank 
r system at equilibrium: 



and consider (rescaled) q r+ i as the new variable. This is allowed only for the classical root 
systems in which r can be any positive integer. Since V of A r+2 has one more element fi r +2 
than that of A r , and As of B r+ \ (BC r+ \) has two more elements e r+ i and — e r +i than that 



W(qi, . . . , q r ) -> W(q r +i) = W(qx, ...,q r , q r +i), 



(6.1) 



of B r (BC r ), we multiply a/ dq r+2 for A r case, 



( V dq r+ i) 2 = dq r+ i for B r (BC r ) case, see 



|£D, Q, O and fl6Tp . 
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6.1 Hermite 

The pre-potential for the A r Calogero system is 

w = -^q 2 + g Yl log( ^' ~ 

l<j<l<r+l 

After rescaling 

V uj 



we obtain from ( 3.1 



~ 1 r+1 fu 

W(z)/g = --z 2 + J^logO - + (z-indep.). 



If we extract a function ip r+ i(z) from 



r+1 I — 

e w/9 ^/dq~^ 2 = (z-indep.) x e" 2 ^ 2 ]J(z - J-qj) x Vd~z 

-i=i V 5 



= (z-indep.) x e~ z ^ 2 H r+1 (z)Vdz 
= i) r+l {z)\fdz, 

it satisfies the orthogonality relation 

dzip n (z)ip m (z) (x 5 nm . 



6.2 Laguerre 

The pre-potential for the B r Calogero system is 

W = --uq 2 + g L Y lo s({Qj ~ Ql){Qj + Qi)j + 9s^2^og 

After rescaling 



1j 

l<j<l<r j=l 



9L 

Qr+1 — \ Z, 

UJ 



we obtain from ( |6.1|) 

W(z)/g L = --z + J^log(z - yJ—Qj) )+-log2;+ (2-indep.), k = 



v 9l 
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2n 



If we extract a function ip r (z) from 

e W ^dq r+1 = (s-indep.) x z^V*/ 2 J]( Z - (^fc)*) * ^ V '^ 
= (z-indep.) x WV 2 / 2 ^^)^ 



= ip r (z)dz, a = k — 1, 
it satisfies the orthogonality relation 

dz ip n (z)4> m (z) oc S n>m . 



6.3 Chebyshev 

This is slightly contrived. The pre-potential for the A r Sutherland system is 

W = g ^ logsin(g i - 

l<j<l<r+l 

By the choice of q ( |5.1| ) and its property q~j = —q r+ 2~j, after defining 

sin q r+2 = z, 

we obtain from ( |6.1| ) 

r+l 

W(z)/g = ^^log(z — sing.,) + (z-indep.). 
If we extract a function ip r+ i(z) from 



r+l 



e w /ff ^/ dq r+2 = (z-indep.) x Y\_( z ~ sing,) x (1 — z 2 ) 1 ^\fdz 

3=1 

= (z-indep.) x (1 - z 2 )- llA T r+l (z)Vd~z 



= ifj r+1 (z)Vdz, 
it satisfies the orthogonality relation 

dzip n (z)ip m (z) oc 5 n . 



38 



6.4 Jacobi 



The pre-potential for the BC r Sutherland system is 

r 

W = g M lo g( sin (<7i - Qi) sin (<?j + qi)j + ^2(gs log sin qj + g L log sin 2%). 

l<j<l<r j=l 

After defining z by 

cos 2q r+ i = z, (6-11) 
we obtain from ( ETfl) (k t = g s /gM, h = 9lI9m) 



W{z)/g M = Y^\og{z - cos 2^) + ^±^l og (l - z) + y log(l + z) + (z-indep.). (6.12) 
If we extract a function tp r {z) from 

r 

e w/9M dq r+1 = (z-indep.) x (1 - z )^+ k ^/ 2 (i + z )>*/* _ cos 2^) x (1 - z 2 )~ 1/2 dz 
= (z-indep.) x (1 - z) a/2 {l + zf' 2 P^\z)dz 

= i[) r (z)dz, a = ki + k 2 — 1, (3 = fc 2 — 1, (6.13) 

it satisfies the orthogonality relation 

dz^ n {z)ijj m {z) oc <J n tn . 
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